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基于 CT 大 度 赋值 的 有 限 元 模型 建立 及 其 在 截 骨 矫 形 中 的 应 用 
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摘要 :目的 探索 一 种 新 的 有 限 元 建 模 方法 以 解决 传统 建 模 方法 在 截 骨 矫形 应 用 中 难以 实现 原 位 CT 灰 度 材 料 赋值 的 问题 。 方 


法 采用 CT 注 层 扫描 方式 采集 1 例 肥 骨 外 旋 扭转 畸形 的 成 人 断层 影像 数据 ,利用 Mimics 软 件 对 采集 的 肥 骨 CT 数据 进行 三 维 重 
线 \ 旋 转角 度 及 参考 平面 并 模拟 截 骨 , 将 稚 骨 后 肥 骨 模型 各 部 导入 有 限 
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行 网 格 划分 并 输出 ,重新 在 Mimics 软 件 中 输入 CT 图 像 的 灰 度 -材料 赋值 关系 公式 ,分 别 对 网 格 模型 进行 非 均 


骨 各 部 分 的 赋值 模型 ,最 后 将 各 部 分 模型 连同 参考 轴 参考 平面 再 次 在 前 处 理 软件 中 完成 截 骨 后 的 装配 组 合 ， 


得 到 截 骨 后 完整 的 肥 骨 有 限 元 模型 ,并 提交 模型 到 有 限 元 求解 器 完成 计算 。 结 果 采用 本 方法 建立 的 肥 骨 截 骨 矫 形 模型 同样 具 


备 基于 CT 灰 度 上 


的 非 均 匀 材 料 属性 ， 


并 可 用 于 


高 效 、 灵 活 和 准确 的 计算 力学 分 析 手 段 。 
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F 截 骨 术 后 结构 的 有 限 元 分 析 。 结 论 本 研究 提出 的 有 限 元 建 模 方法 既 保留 了 基于 
CT 灰 度 的 材料 赋值 方法 的 准确 性 ,又 能 解决 传统 赋值 方法 难以 实现 的 骨骼 模型 部 件 之 间 空 间 变换 问题 ,从 而 为 矫形 外 科 提 供 


Finite element modeling of material property assignment based on CT gray value and its 


application in simulation of osteotomy for deformities 
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‘Department of Anatomy, School of Basic Medical Sciences, Southern Medical University, Guangzhou 510515, China; “Department of 
Orthopedic Surgery, Third Affiliated Hospital of Southern Medical University, Guangzhou, 510630, China 


Abstract: Objective To explore a new method for finite element modeling to achieve material property assignment based on in 
situ CT gray value in simulated osteotomies for deformities. Methods A CT scan dataset of the lower limb of a patient with 
extorsion deformity was obtained for three-dimensional reconstruction using Mimics software and preparing a solid model. mm 
the CAD software, the parameters for osteotomy simulation were defined including the navigation axis, rotation angle and 
reference plane. The tibia model was imported to the FEA pre-processing software for meshing procedure and then exported 
to Mimics. All the segments of the tibia meshed model were assigned uneven material Properties based on the relationship 
between CT gray values and material properties in the Mimics software. Finally, all the segments of the tibia model, reference 
axis and reference plane were assembled in the pre-processing software to form a full finite element model of a corrected tibia, 
which was submitted to resolver for biomechanical analysis. Results The tibia model established using our modeling method 
had inhomogeneous material properties based on CT gray values, and was available for finite element analysis for the 
simulation of osteotomy. Conclusions The proposed finite element modeling method, which retains the accuracy of the 
material property assignment based on CT gray value, can solve the reposition problem commonly seen in modeling via the 
routine method of property assignment and provides an efficient, flexible and accurate computational biomechanical analysis 
method for orthopedic surgery:. 
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随 着 数字 医学 技术 与 计算 机 硬件 技术 的 快速 发 展 ， 
有 限 元 法 作为 一 种 离散 数值 计算 方法 已 被 广泛 用 于 骨 
科 生物 力学 领域 "。 在 过 去 数 十 年 中 ,如 何 获得 更 为 合 


理 的 材料 属性 赋值 和 更 为 精确 的 几何 模型 一 直 是 骨科 
生物 力学 有 限 元 分 析 的 热点 问题 ”。 目 前 ,对 于 骨骼 有 
限 元 建 模 的 材料 参数 赋值 方法 主要 有 两 种 ,一 种 是 采用 
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复合 结构 建 模 的 均一 赋值 方法 , 另 一 种 是 基于 CT 灰 度 
经 验 公 式 的 非 均匀 赋值 法 “。 前 者 将 皮质 骨 和 松 质 骨 
分 别 简化 作 均 一 属性 的 材料 ,赋值 过 程 简便 ,快速 ,可 对 
模型 进行 各 种 空间 位 置 变换 操作 ,但 模型 的 简化 赋值 可 
能 导致 最 终 分 析 结 果 准确 性 不 足 ; 而 后 者 属于 原 位 赋值 
方法 ,具有 个 体 化 建 模 的 特点 ,但 难以 进行 模型 的 虚拟 
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操作 (如 截 骨 矫形 \ 虚 拟 复位 等 )”。 然 而 ,同时 兼 具 上 述 
两 种 方法 优点 的 有 限 元 建 模 方 法 在 国内 未 见报 道 。 为 
此 ,本 文 以 肥 骨 外 旋 扭 转 畸 形 病例 数据 为 应 用 范例 , 探 
索 一 种 新 的 有 限 元 建 模 方 法 以 满足 截 骨 矫形 仿真 分 析 
的 需求 。 


1 材料 和 方法 
1.1 数据 采集 及 肥 骨 实体 模型 的 三 维 重 建 

选取 1 名 成 人 石 侧 肥 骨 外 旋 扭 转 畸 形 病 例 作为 数 
据 采 集 来 源 ,采用 64 排 螺旋 CT 机 (西门 子 , 德 国 ) 沿 肥 骨 
于 长 轴 方 向 进行 CT 扫描 ,范围 覆盖 膝 关 节 及 躁 关节 , 扫 
描 层 距 为 1 mm, 共 获得 连续 横断 面 1000 张 5312x512 像 
素 的 CT 图 像 ,并 将 获得 的 CT 图 像 以 通用 的 DICOM 标 
准 格式 导 和 人 Mimics 14.0 软件 (Materialise, 比利时 ) 进 
行 肥 骨 模型 三 维 重 建 ,参考 既往 文献 的 标准 测量 方 
法 在 三 维 模型 上 测 得 肥 骨 外 旋 角 度 为 41.8° ,符合 所 
转 畸 形 诊断 的 解剖 参数 标准 ,导出 为 STL 文 件 保存 。 
将 STL 文 件 导 入 逆向 工程 软件 Geomagic Studio 2012 
(Geomasgic, 美 国 ) 中 进行 降 噪 和 光 顺 表面 处 理 , 进 而 构 
建 出 肥 骨 实体 模型 ,过 程 如 图 1 所 示 。 


Jsurface modification 
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图 1 肥 骨 三 维 模型 的 表面 处 理 与 实体 建 模 
Fig.1 Surface treatment and solid modeling of the 3D tibia 
model. 


1.2 盈 骨 模型 虚拟 截 骨 矫 形 及 参考 轴 面 的 定义 

将 肥 骨 实体 模型 导入 三 维 建 模 软 件 UG NX 8.5( 西 
门 子 , 德 国 ) 中 确定 肥 骨 上 部 干 骨 端的 截 骨 平面 并 以 此 
建立 参考 平面 ,同时 以 拟 合 的 肥 骨 长 轴 作 为 截 骨 旋 转轴 
并 建立 相应 的 空间 直线 ,以 参考 平面 对 肥 骨 进行 虚拟 截 
骨 操作 ,获得 截 骨 后 上 下 两 个 部 件 子 模型 ,并 以 此 作为 
划分 单元 网 格 的 对 象 (图 2)。 
1.3 CT 灰 度 赋值 法 有 限 元 模型 建立 

将 旋转 操作 前 的 各 部 肥 骨 子 模型 导入 有 限 元 分 析 
软件 Abaqus 6.14( 达 索 , 法 国 ) 中 划分 网 格 后 并 以 inp 文 
件 格 式 导 入 Mimics 软件 中 进行 材料 属性 赋值 , 设 定 图 
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图 2 肥 骨 截 骨 平面 的 建立 (4) 与 虚拟 截 骨 操 作 (B) 
Fig.2 Simulation of osteotomy on the tibia model. A: 
Definition of the plane of osteotomy; B: Process of 
simulative osteotomy. 


像 CT 值 HU 与 表 观 密度 、 表 观 密度 与 弹性 模 量 的 关系 
公式 “后 ,完成 基于 CT 图 像 灰 度 的 原 位 赋值 过 程 ,最 后 
以 inp 格 式 把 肥 骨 上 下 部 分 子 模型 导 回 至 有 限 元 软件 前 
处 理 模块 中 (图 3)。 


> 3 
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3 基于 CT 灰 度 -材料 属性 关系 公式 赋予 脏 骨 非 均 匀 材 料 属性 


Fig.3 Uneven material property assignment of the tibia based 


on the relationship between CT gray values and material 
properties. 


1.4 肥 骨 虚拟 旋转 截 骨 位 置 变 换 与 有 限 元 前 处 理 

将 先前 定义 好 的 参考 轴 模 型 导入 Abaqus 软件 中 ， 
在 装配 模式 下 按照 预先 规划 的 10°" 内 旋 角 度 , 对 肥 骨 下 
部 模型 沿 着 参考 旋转 轴 进 行 旋转 位 置 变换 操作 ,对 截 骨 
平面 设置 为 Tie 约 束 ,得 到 调整 位 置 后 完整 的 肥 骨 旋转 
截 骨 重 装配 模型 (图 4)。 参 照 文献 数据 ,在 肥 骨 上 端 施 
加 20 Nm 的 招 算 载荷 ,固定 肥 骨 远 端 躁 关节 平面 作为 模 
型 边界 条 件 , 最 后 提交 至 求解 右 计 算 分 析 截 骨 后 的 应 力 和 
位 移 分 布 结 果 。 至 此 ,整个 建 模 分 析 过 程 全 部 完成 (图 5)。 


2 结果 

通过 CT 灰 度 值 赋 值 法 建立 的 肥 骨 旋转 截 骨 矫 形 
有 限 元 模型 可 顺利 提交 分 析 , 并 得 到 肥 骨 旋转 截 骨 后 的 
应 力 和 位 移 分 布 结果 (图 6)。 在 用 骨 上 有 段 施加 20 Nm 
的 扭矩 载荷 作用 下 ,对 于 肥 骨 的 应 力 分 布 方面 ,旋转 截 
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图 4 完成 材料 属性 赋值 后 的 用 骨 虚拟 旋转 截 骨 操 作 
Fig.4 Simulated rotational osteotomy of the tibia after 
properties assignment. 
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图 5 建 模 分 析 方 法 流程 图 
Fig.5 Flow chart of modeling and analysis. 


结果 可 靠 性 的 保证 。 本 人 研究 采用 了 非 均 质 材料 赋值 的 
建 模 方法 构建 了 肥 骨 有 限 元 模型 ,在 理论 上 更 为 接近 人 
体 骨骼 材料 的 力学 特性 ,可 使 模型 的 分 析 计算 具有 较 高 
的 可 靠 性 。 这 一 方法 在 国内 外 的 文献 报道 中 均 已 得 到 
了 应 用 ,同时 也 得 到 了 体外 生物 力学 实验 的 文献 数据 支 
持 ”"。 相 比 之 下 ,以往 也 有 部 分 人 体 骨 骼 的 有 限 元 研 
究 对 骨骼 的 材料 属性 作 均 匀 分 布 简化 处 理 , 将 骨骼 模型 
简单 划分 为 密 质 骨 和 松 质 骨 两 种 材料 进行 赋值 "。 该 
方法 仅 包 括 了 三 维 重建 .虚拟 装配 和 有 限 元 前 处 理 3 个 
部 分 ,在 一 定 程 度 上 提高 了 分 析 建 模 的 效率 。 然 而 ,对 
于 骨 密 度 差 异 较 小 的 模糊 交界 区 域 采 用 均一 赋值 处 理 
将 在 一 定 程 度 上 人 为 提高 或 降低 整体 骨髓 的 结构 刚度 ， 
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骨 后 的 应 力 分 布 云图 显示 (图 6A 工 ~ 工 ) ,最 大 应 力 峰 值 
主要 分 布 在 肥 骨 旋转 截 骨 平面 部 位 ,应 力 峰 值 为 306.8 
MPa; 当 忽 略 截 骨 平 面部 位 后 (图 6A 卫 ~WV ), 肥 骨 最 大 
应 力 分 布 在 肥 骨 中 下 段 , 最 大 应 力 峰值 为 56.83 MPa; 
对 于 位 移 分 布 方面 (图 6B I ~JV ) ,同样 扭矩 作用 载荷 
下 , 肥 骨 的 最 大 位 移 分 布 主要 位 于 肥 骨 上 有 段 , 最 大 合 位 
移 峰 值 为 2.15 mm, 其 中 X 和 YY 方向 位 移 分 量 为 主要 位 
移 分 量 , 最 大 峰值 分 别 为 1.063 mm 和 2.005 mm,Z 方 向 
位 移 分 量 的 最 大 峰值 为 0.123 mm。 


3 讨论 

在 建 模 方法 方面 ,有 限 元 分 析 的 准确 性 很 大 程度 上 
取决 于 建 模 过 程 的 合理 性 ,而 骨骼 材料 属性 的 赋值 是 骨 
骼 有 限 元 分 析 中 的 关键 环节 ,也 是 有 限 元 模型 计算 分 析 
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从 而 影响 最 终 的 分 析 结 果 精 度 ,对 比 骨 骼 真实 应 变 情况 
容易 产生 较 大 偏差 ,这 将 给 不 同 模型 分 组 之 间 的 有 限 元 
对 比分 析 结 果 带 来 了 显著 的 不 稳定 性 "; 而 本 人 研究 采用 
的 非 均匀 材料 赋值 方法 涉及 了 多 个 建 模 软件 之 间 的 文 
件 交 互 对 接 , 步 又 程序 相对 细 化 ,尽管 影响 了 建 模 过 程 
的 效率 ,但 可 获得 更 为 接近 人 体 骨髓 真实 材料 性 质 的 模 
型 ,有 效 提高 有 限 元 分 析 结 果 的 稳定 性 。 

在 Mimics 软 件 中 ,用 户 能 根据 骨骼 CT 图 像 中 不 同 
区 域 的 灰 度 值 按 经 验 公式 赋予 模型 非 均 匀 的 材料 属 
性 ”5。 非 均匀 赋值 方法 可 以 全 面 分 析 骨 组 织 复杂 的 
密度 变化 情况 ,基于 此 法 的 有 限 元 分 析 结 果 具 有 更 高 的 
准确 度 "“。 然 而 ,这 一 方法 要 求 网 格 模型 空间 位 置 须 
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图 6 扭矩 作用 下 胖 骨 应 力 分 布 云 图 (4 ) 和 胖 骨 位 移 分 布 云 图 (B) 


Fig.6 Results of finite element analysis of tibia under torque. A: Distribution of stress; B: Distribution of displacement. 


和 原 有 CT 图 像 覆 盖 范 围 保 持 一 致 , 当 模型 空间 位 置 发 
生变 化 时 , 原 位 CT 灰 度 赋予 材料 属性 的 方法 将 完全 失 
效 ,因此 应 用 既往 方法 将 无 法 在 截 骨 矫形 后 的 骨髓 有限 
元 分 析 中 应 用 。 因 此 ,本 研究 提出 的 建 模 技术 路 线 , 有 
效 解决 了 既往 方法 不 能 直接 对 模型 位 置 变换 后 的 赋值 
问题 。 笔 者 以 肥 骨 扭转 畸形 截 骨 分 析 为 应 用 实例 ,利用 
不 同 软件 之 间 的 模型 文件 转换 接口 ,实现 了 材料 原 位 赋 
值 和 虚拟 截 骨 之 间 的 兼容 ,形成 一 套 新 的 有 限 元 建 模 流 
程 。 尽 管 国 外 已 有 研究 机 构 自 主 研发 了 相应 的 映射 配 
对 方法 来 解决 这 一 问题 ,该 方法 尚 不 能 很 好 地 解决 局 部 
网 格 单元 缺失 的 问题 "。 因 此 ,本 研究 的 技术 方案 相对 
更 为 简单 高 效 ,适用 性 较 强 。 

在 分 析 结 果 方 面 ,本 研究 方法 构建 的 有 限 元 模型 可 
顺利 提交 分 析 , 从 应 力 和 位 移 分 布 的 云图 来 看 ,计算 结 
果 较 为 接近 临床 实际 情况 ,扭转 刚度 数值 也 较为 接近 以 
往 文献 数据 ”。 从 整个 建 模 过 程 上 看 ,本 研究 基本 实现 
了 上 肥 骨 旋转 截 骨 矫形 术 后 的 生物 力学 预测 分 析 , 在 保证 


各 部 份 骨骼 按 CT 灰 度 值 赋值 法 进行 非 均匀 材料 属性 赋 
值 的 前 提 下 , 截 骨 后 的 骨骼 各 部 也 能 根据 术 前 设计 参考 
轴 面 及 角度 实现 准确 的 旋转 截 骨 , 从 而 满足 了 骨骼 矫形 
手术 的 有 限 元 仿真 分 析 需 求 。 从 结果 上 看 ,本 研究 在 对 
有 径 骨 上 段 施加 扭矩 载荷 的 作用 下 ,在 截 骨 平面 部 位 产生 
了 较 大 的 应 力 集中 现象 ,最 大 峰值 达到 了 306.8 MPa， 
造成 这 一 结果 的 原因 主要 是 因为 本 研究 并 未 考虑 截 骨 
术 后 早期 的 内 固定 器 械 对 模型 支撑 效应 ,同时 未 对 规 骨 
矫形 的 后 期 截 骨 平面 骨 愈合 形态 进行 模拟 ,因此 仅仅 对 
截 骨 平面 的 相 邻 单元 节点 进行 了 绑 定 约束 ,从 而 形成 了 
截 骨 平面 单元 节点 之 间 形 成 了 刚性 连接 ,使 得 结果 出 现 
了 较 高 的 应 力 峰值 。 通 过 忽略 截 骨 面 高 应 力 单元 可 消 
除 截 骨 面 绑 定 约束 设 定 对 结果 的 影响 ,由 此 获得 的 应 力 
分 布 特点 更 为 合理 ,而 大 体 应 力 数值 水 平 也 与 文献 报道 
同等 载荷 下 的 分 析 结 果 基 本 一 臻 。 因 此 ,本 模型 分 析 结 
果 基 本 符合 临床 实际 情况 。 而 在 实际 分 析 中 ,应 用 本 方 
法 建 模 时 应 明确 模型 是 处 于 截 骨 矫形 术 后 的 哪 一 个 时 
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期 ,明确 是 否 存 在 相应 的 植 和 信物 共同 承担 生理 载荷 和 和 骨 
愈合 情况 ,并 根据 实际 情况 调整 模型 设 定 。 

总 的 来 说 ,本 研究 提出 的 建 模 方法 在 应 用 拓展 性 方 
面 显示 出 一 定 的 优势 :首先 ,本 方法 可 用 于 不 同 骨 密度 
水 平 的 活体 个 体 化 矫形 外 科 生 物力 学 预测 分 析 ; 其 次 ， 
实施 过 程 中 采用 的 建 模 软件 均 为 现 有 可 用 的 商业 建 模 
软件 ,不 需要 自行 开发 空间 映射 配对 赋值 软件 ,简化 了 
建 模 过 程 ,具有 较 高 的 可 操作 性 。 然 而 ,本 研究 在 模型 
的 定量 评价 方面 也 存在 着 一 定 的 局 限 性 。 目 前 ,模型 的 
定量 评价 主要 有 同期 生物 力学 实验 验证 文献 数据 验证 
和 临床 验证 3 种 方式 。 第 1 种 方式 是 利用 尸体 标本 的 实 
验 生 物力 学 数据 对 模型 进行 定量 评估 ,而 在 本 研究 中 建 
芯 的 肥 骨 模型 源 于 扭转 畸形 患者 的 CT 扫描 数据 ,因此 
生物 力学 实验 无 法 在 该 患者 肥 骨 上 实施 。 而 获取 骨髓 
畸形 的 尸体 标本 极为 困难 ,同样 导致 实验 定量 评估 难以 
开展 。 对 于 第 2 种 方式 ,与 本 人 研究 相关 的 肥 骨 扭转 畸形 
分 析 鲜 有 文献 报道 ,上 且 既 往 研 究 中 的 个 体 和 本 研究 中 的 
患者 个 体 存在 着 一 定 的 解剖 差异 ,因而 可 比 性 受到 一 定 
的 影响 。 而 本 研究 选取 了 同样 载荷 条 件 下 的 正常 胜 骨 
分 析 文 献 数 据 进行 对 比 ,尽管 正常 肥 骨 与 扭转 畸形 肥 骨 
存在 着 结构 上 的 差异 ,但 本 研究 中 整体 肥 骨 的 应 力 分 布 
特点 和 数值 基本 处 于 合理 范围 ,这 在 一 定 程度 上 验证 了 
模型 的 有 效 性 。 而 对 于 第 3 种 方式 ,本 肥 骨 模型 的 高 应 力 
区 主要 位 于 胜 骨 截面 发 生 明 显 变 化 的 中 下 1/3 部 位 ,该 
部 位 是 临床 上 胜 骨 骨折 发 生 的 最 常见 部 位 ,说 明 模型 大 
体 符合 临床 实际 情况 。 而 在 后 续 的 研究 中 ,采用 活体 动 
物 骨骼 畸形 造 模 的 方法 有 望 解决 上 述 定量 评估 问题 ,但 
活体 验证 实施 时 间 较 长 , 且 受 物种 之 间 解 剖 结构 差异 的 
影响 ,同样 存在 一 定 的 局 限 性 。 而 通过 活体 的 放射 立体 
照相 测量 分 析 (RSA) 和 有 限 元 分 析 相 结合 的 方式 ,有 望 
为 模型 临床 评估 提供 一 种 更 为 理想 的 手段 ”。 

本 研究 提出 了 一 种 基于 现 有 医学 三 维 建 模 软件 的 
骨骼 有 限 元 建 模 方法 ,并 以 一 例 肥 骨 扭转 畸形 病例 的 旋 
转 截 骨 矫 形 有 限 元 分 析 介 绍 了 该 方法 的 应 用 。 结 果 表 
明 , 该 方法 实现 了 骨髓 CT 图 像 灰 度 - 材 料 参数 赋值 的 同 
时 ,还 可 进行 模型 空间 位 置 变换 操作 ,因而 可 为 临床 截 
骨 矫 形 手 术 的 有 限 元 分 析 预 测 提供 简单 .可靠 和 准确 的 
建 模 手段 。 
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